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We point out that for solar neutrino oscillations with the mass-squared difference of Am 2 



10" 



10 eV , traditionally known as "vacuum oscillation" range, the solar matter effects are 



non-negligible, particularly for the low energy pp neutrinos. One consequence of this is that the 
values of the mixing angle 9 and n/2 — 9 are not equivalent, leading to the need to consider the 
entire physical range of the mixing angle < 9 < n/2 when determining the allowed values of the 
neutrino oscillation parameters. 



1. The field of solar neutrino physics is currently un- 
dergoing a remarkable change. For 30 years the goal 
was simply to confirm the deficit of solar neutrinos. The 
latest experiments, however, such as Super-Kamiokande, 
SNO, Borexino, KamLAND, etc, aim to accomplish more 
than that. By collecting high statistics real-time data 
sets on different components of the solar neutrino spec- 
trum, they hope to obtain unequivocal proof of neu- 
trino oscillations and measure the oscillation parameters. 
With the physics of solar neutrinos quickly becoming a 
precision science, it is more important then ever to ensure 
that all relevant physical effects are taken into account 
and the right parameter set is used. 

It has been a long-standing tradition in solar neutrino 
physics to present experimental results in the Am 2 — 
sin 2 29 space and to treat separately the "vacuum os- 
cillation" (Am 2 - 1CT 11 - 1CT 9 eV 2 ) and the MSW 
(Am 2 ~ 10~ 8 — 10~ 3 eV 2 ) regions. In the vacuum oscil- 
lation region the neutrino survival probability (i.e. the 
probability to be detected as v e ) was always computed 
according to the canonical formula, 



P = 1 - sin 2 29 sin 2 1.27 
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E 



(1) 



where the neutrino energy E is in GeV, the distance L in 
km, and the mass-squared splitting Am 2 in eV 2 . Eq. (|l|) 
makes sin 2 29 seem like a natural parameter choice. As 
sin 2 29 runs from to 1, the corresponding range of the 
mixing angle isO < 9 < ir/4. There is no need to 
treat separately the case of Am 2 < (or equivalently 
7r/4 < 9 < 7r/2), since Eq. (|l|) is invariant with respect 
to Am 2 -> -Am 2 (9 -> tt/2 - 9). 

The situation is different in the MSW region, since 
neutrino interactions with matter are manifestly flavor- 
dependent. It is well known that for |Am 2 | > 1CP 8 eV 2 
matter effects in the Sun and Earth can be quite large. 
In this case, if one still chooses to limit the range of 
the mixing angle to < 9 < 7r/4, one must consider 
both signs of Am 2 to describe all physically inequivalent 
situations. As was argued in M, to exhibit the continuity 
of physics around the maximal mixing, it is more natural 
to keep the same sign of Am 2 and to vary the mixing 
angle in the range < 9 < w/2. 



Historically, a possible argument in favor of not consid- 
ering 9 > 7r/4 in the MSW region might have been that 
this half of the parameter space is "uninteresting" , since 
for 9 > 7r/4 there is no level-crossing in the Sun and the 
neutrino survival probability is always greater than 1/2. 
However, a detailed analysis reveals that allowed MSW 
regions can extend to maximal mixing and beyond, as 
was explored in j| (see also || and Q for a treatment of 
3- and 4- neutrino mixing schemes). 

In this letter we point out that for solar neutrinos with 
low energies, particularly the pp neutrinos, the solar mat- 
ter effects can be relevant even for neutrino oscillations 



with Ar 



10 



-10 



10" 9 eV 2 . These effects break the 



symmetry between 9 and n/2 — 9 making it necessary 
to consider the full physical range of the mixing angle 
< 9 < tt/2 even in the "vacuum oscillation" case. 

2. For simplicity, we will only consider here the two- 
generation mixing. If neutrino masses are nonzero then, 
in general, the mass eigenstates 1 1^1,2) a- re different from 
the flavor eigenstates \v e .fi)- The relationship between 
the two bases is given in terms of the mixing angle 9: 



\v\) = cos9\v e ) — sin#|z^}, 
1 1/2) = sn\9\v e ) + cos0|i/ M ). 



(2) 



In our convention [1/2) is always the heavier of the two 
eigenstates, i.e. Am 2 = m\ — m\ > 0. Then, as already 
mentioned, < 9 < tt/2 encompasses all physically dif- 
ferent situations. 

Neutrinos are created in the Sun's core and exit the 
Sun in the superposition of \v±) and 1 1^2 ) - For Am 2 in 
the vacuum oscillation region, the neutrino is produced 
almost completely in the heavy Hamiltonian eigenstate 
\v+). In this case, if the evolution inside the Sun is adia- 
batic, the exit state is purely \v2)- In the case of a nona- 
diabatic transition there is also a nonzero probability P c 
to find the neutrino in the \v\) state (a "level crossing" 
probability). For a given value of P c , the survival prob- 
ability for neutrinos arriving at the Earth is determined 
by simple 2-state quantum mechanics [ 



P = P. cos 2 



(1 - P r ) sin 2 



2y/Pc{\ - P c ) sin 9 cos 9 cos 2.54 
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Here S is a phase acquired when neutrinos traverse the 
Sun. In our analysis it is determined numerically 
Units are the same as in Eq. (|j). 

In the adiabatic limit P c = and Eq. (|J) yields 
P = sin 2 9. Neutrinos exit the Sun in the heavy mass 
eigenstate and do not oscillate in vacuum. In the op- 
posite limit of small Am 2 , when the neutrino evolution 
in the Sun is "extremely nonadiabatic" , P c — > cos 2 9. It 
is trivial to verify that Eq. (||) in this limit reduces to 
Eq. (Q). It has been assumed that in the vacuum oscil- 
lation region this limit is reached. Remarkably, however, 
this is not always the case for the low energy solar neu- 
trinos, especially the pp neutrinos (E v < 0.42 MeV). 

The most reliable way to compute P c is by numerically 
solving the Schrodinger equation in the Sun for different 
values of Am 2 and 9. We do this using the latest avail- 
able BP2000 solar profile |5). Fig. [I] shows contours of 
constant P c for the energy of 7 Be neutrino (solid lines). 
Note that the variable on the horizontal axis is tan 2 9. 
With this choice, points 9 and ir/2 — 9 are located sym- 
metrically on the logarithmic scale about tan 2 9=1 (see 
pof) The figure demonstrates that the contours are 
not symmetric with respect to the tan 2 9 = 1 line, ex- 
cept in the region of Am 2 /E u < 10~ 10 eV 2 /MeV, where 
the extreme nonadiabatic limit is reached. This simple 
observation is the crucial point of this letter. 
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FIG. 1. Contours of constant level crossing probability P c 
for neutrino energy of 0.863 MeV ( 7 Be line). The solid lines 
are the results of numerical calculations using the BP2000 
solar profile. The dashed lines correspond to using the expo- 
nential profile formula with ro = -Rq/18.4 = 3.77 x 10 4 km. 

In the MSW region the value of P c is often computed 
using the analytical result p| 



Pr 



e 7 cos 2 e _ i 
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(4) 
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valid for the exponential solar profile n e cx exp(— r/ro), 
with r = i?©/10.54 = 6.60 x 10 4 km Q. Although orig- 
inally derived for 9 < tt/4, it also applies when 9 > 7r/4, 
as was demonstrated in jlj . In the region relevant for vac- 
uum oscillation, however, 0.9i?o < R < R@, the profile 
falls off faster than the exponential with ro = i?0/lO.54. 
Nevertheless, Eq. ([|) can still be used with the appro- 
priately chosen value of ro. The dashed lines in Fig. [l] 
show the contours of P c computed using Eq. (jj) with 
r = Rq/18A = 3.77 x I0 4 km. As can be seen from the 
figure, the agreement between the two sets of contours 
for Am 2 <4x 10~ 9 eV 2 is very good. Note that a sim- 
ilar result was arrived at in juj for 9 < tt/4, where the 
value of ro = R Q x 0.065 = 6.5 x 10 4 km was obtained. 

Fig. ^ can also be used to read off the values of P c 
for different neutrino energies, since P c depends on E v 
through the combination Am 2 /E„. It is obvious that for 
neutrinos of lower energies P c starts deviating from its 
"extreme nonadiabatic" value at even smaller values of 
Am 2 , and vice versa. Consequently, as will be seen later, 
the solar matter effects on vacuum oscillations are most 
important at the gallium experiments, which are sensi- 
tive to the pp neutrinos, while the Super-Kamiokande 
experiment is practically unaffected. 

Using Eqs. (|4|,||), it is possible to derive a corrected 
form of Eq. (|Tj) , by retaining in the expansion terms linear 
in 7: 



P= 1 - 

+ 0( 7 2 



7 \ 9 o ( Am L 
1 + 1 cos 29) sin 2 29 sin 2 ( 1.27—— 
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where 



Notice that the first order correction contains cos 29 
and hence is manifestly not invariant under the transfor- 
mation 9 -> n/2-6. Using Eq. (|) withr = 3.8xl0 4 km, 
we see that for the pp neutrinos [E v < 0.42 MeV) this 
correction is indeed non-negligible already for Am 2 ~ 
10- 10 - 10- 9 eV 2 . 

With matter effects being relevant already at Am 2 > 
10~ 10 eV 2 one might wonder if the separation between 
vacuum oscillation solutions and MSW solutions is some- 
what artificial. To fix the terminology, we will adopt a 
definition of vacuum oscillations as the situation when 
the value of neutrino survival probability depends on the 
distance L from the Sun, regardless of whether matter 
effects are negligible or not. The transition between the 
vacuum and the MSW regions will be discussed shortly. 

3. To illustrate the role of matter effects in vacuum 
oscillations, we present fits to the total rates of the Home- 
stake H, GALLEX H| and SAGE @, and Super- 
Kamiokande fl§| l experiments. We combine experimental 
rates and uncertainties for the two gallium experiments 
and use the latest available 825-day Super-Kamiokande 
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FIG. 2. Regions allowed by total rates of GALLEX and SAGE only (A), GALLEX, SAGE, and Super-Kamiokande (B), and 
GALLEX, SAGE, Homestake, and Super-Kamiokande (C). Black outlines correspond to neglecting the solar matter effects. 



data set. The experimental results are conveniently col- 
lected and tabulated in ||] . 

We fit the data to the theoretical predictions of the 
BP98 standard solar model [pTsf] . Predicted fluxes and 
uncertainties for various solar reactions were kindly made 
available by J. N. Bahcall at |2(J. To compute the rate 
suppression caused by neutrino oscillations, we numeri- 
cally integrate the neutrino survival probability, Eq. (|^), 
over the energy spectra of the pp, 7 Be, 8 B, pep, 13 N, and 
15 neutrinos. In addition, to account for the fact that 
the Earth-Sun distance L varies throughout the year as 
a consequence of the eccentricity of the Earth's orbit 

L = Lq(1 — e cos(2-7rf/year)) (7) 

we also integrate over time to find an average event rate. 
In Eq. (0) t is time measured in years from the perihelion, 
Lq = 1.5 xlO 8 km is one astronomical unit, and e = 1.7%. 

In Fig. ^| (A) we show the vacuum oscillation regions 
allowed by the total rates of GALLEX and SAGE. For 
comparison, we also show the regions one would obtain by 
neglecting the neutrino interactions with the solar mat- 
ter (dark outlines), i.e. by setting P c = cos 2 8 (black 
contours). The allowed regions were defined as the sets 
of points where the theoretically predicted and experi- 
mentally observed rates are consistent with each other 
at the 2cr CL. for 1 d.o.f. (x 2 = 4.0) The plot 

demonstrates that the matter effects at the gallium ex- 
periments are quite important, with their contribution 
being significant for Am 2 >2x 10~ 10 eV 2 . 

The remaining two plots in Fig. || show the vacuum 
regions allowed at 3cr CL. by the rates of GALLEX, 
SAGE, and Super-Kamiokande (B) (2 d.o.f., X 2 = H-83, 
in the same convention as before), and all four experi- 
ments combined (C) (3 d.o.f., \ 2 — 14.15). In order to 
properly account for the correlation between the theoret- 
ical errors of the different experiments, we followed the 
technique developed in |2^] and ||. The matter effects 
are noticeable for Am 2 > 6 x 10~ eV 2 . 
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FIG. 3. The sensitivity region of the Borexino experiment 
to anomalous seasonal variations for the full range of the mix- 
ing angle (95% C.L.). Notice the asymmetry for large Am 2 . 

4. An important question is how well future experi- 
ments will be able to cover vacuum oscillation solutions 
with 9 > 7r/4. In Fig. || we show the sensitivity of 
the Borexino experiment to anomalous seasonal varia- 
tions for the entire physical range of the mixing angle 
< < tt/2. This is an extension of the analysis per- 
formed in Q , where the details of the procedure are de- 
scribed. The sensitivity region shows a clear asymmetry 
as a result of the solar matter effects. 

Fig. [3] gives us an opportunity to discuss the extent 
of the vacuum oscillation region. There are two primary 
physical reasons why the neutrino event rate becomes in- 
dependent of L (and anomalous seasonal variations dis- 
appear) for sufficiently large Am 2 : 

• Adiabatic evolution in the Sun. As P c — > the last 
term in Eq. (||) vanishes. 
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• Integration over neutrino energy spectrum. To 
compute the event rate one has to integrate Eq. (||) 
over neutrino energies. For sufficiently large Am 2 
the last term averages out to zero, leading effec- 
tively to the loss of coherence between the two mass 
eigenstates. 

As Am 2 increases, coherence is first lost for reac- 
tions with broad energy spectra, such as pp and 8 B, 
and persist the longest for neutrinos produced in two- 
body final states. The most important such reaction is 
7 Be+e~ — > 7 Li +v e , which produces the 7 Be neutrinos. 
The 7 Be neutrinos have an energy spread of only a few 
keV, arising from the Doppler shift due to the motion 
of the 7 Be nucleus and the thermal kinetic energy of the 
electron. A detailed discussion of this phenomenon can 
be found in @,|. 

In order to properly take these effects into account, in 
our codes we numerically integrate over the exact 7 Be 
line profile, computed in ^3j. As Fig. || shows, the neu- 
trino survival probability becomes independent of L for 
Am 2 > 6 x 10 -9 eV 2 . For this reason, we present our 
fits for Am 2 ranging from 10 _u cV 2 to 10~ 8 eV 2 . Un- 
fortunately, in the literature vacuum oscillations are usu- 
ally studied in the range from 10 -11 eV 2 to 10~ 9 eV 2 
pi] , p5p |], although the allowed regions in all these pa- 
pers seem to extend above 10 -9 eV 2 . 

5. In summary, the preceding examples clearly illus- 
trate the importance of including the solar matter ef- 
fects when studying vacuum oscillation of solar neutri- 
nos with Am 2 > 10 -10 eV 2 . Because to describe such 
effects one has to use the full range of the mixing angle 
< 6 < 7r/2, future fits to the data should be extended 
to 9 > 7r/4. This seems especially important in light 
of the latest analyses |2f|, jij], which in addition to the 
total rates also use the information on the neutrino spec- 
trum and time variations at Super-Kamiokande. In this 
case the allowed vacuum oscillation regions are mostly 
located in the Am 2 > 4 x 10~ 10 cV 2 region {|, pre- 
cisely where the matter effects are relevant. (The best 
fit to the Super-Kamiokande electron recoil spectrum is 
achieved for Am 2 = 6.3 x KT 10 eV 2 , sin 2 20 = 1 ||.) 
It would be very desirable to repeat these analyses with 
the solar matter effects included. 

Additionally, since the 7 Be neutrinos remain (par- 
tially) coherent for Am 2 > 1CU 9 cV 2 , it is desirable to 
present the results of the fits in the range fCT 11 eV 2 < 
Am 2 < fCT 8 eV 2 , as was done in 101. 
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